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Suppressing Proton Decay in the Minimal SO(IO) Model 



O 

o 

(N 
o 

CD 

Q 

00 



> 
o 

O 



Oh: 

CD 



X 



Bhaskar Dutta*, Yukihiro Mimura*, and R.N. Mohapatra''' 

"Department of Physics, University of Regina, Regina, Saskatchewan S4S OAS, Canada 
''Department of Physics, University of Maryland, College Park, MD 20742, USA 

(Dated: February 2, 2008) 

We show that in a class of minimal supersymmetric SO(IO) models which have been found to 
be quite successful in predicting neutrino mixings, all proton decay modes can be suppressed by a 
particular choice of Yukawa textures. This suppression works for contributions from both left and 
right operators for nucleon decay and for arbitrary tan j3. The required texture not only fits all 
lepton and quark masses as well as CKM parameters but it also predicts neutrino mixing parameter 
Ue3 and Dirac CP phase sin|(5A/jvs| to be 0.07-0.09 and 0.3-0.7 respectively. We also discuss the 
relation between the GUT symmetry breaking parameters for the origin of these textures. 



INTRODUCTION 

Recent experiments in neutrino physics which have es- 
tablished that neutrino masses are in the sub-eV range 
have provided a new reason for taking supersymmetric 
grand unified theories (SUSY GUTs) seriously. This has 
to do with using the seesaw mechanism to understand 
the extreme smallness of neutrino masses compared to 
charged fermion masses [ij. The seesaw mechanism con- 
sists of extending the standard model by introducing 
three right-handed neutrinos with large Majorana masses 
Mfl. Simple arguments based on the atmospheric oscil- 
lation data then tell us that at least one Mji is close to 
the conventional GUT scale AIjj ^ 10^® GeV, perhaps 
suggesting that the seesaw scale, Mb, and GUT scale are 
one and the same. 

The minimal GUT that unifies the right-handed neu- 
trinos with the other fermions of the standard model and 
leads to the seesaw mechanism is the SO(IO) model where 
all quarks and leptons are unified into one 16 dimensional 
spinor multiplet. This raises the possibility that the va- 
riety of masses and mixings of quarks and leptons can 
be understood in terms of a smaller number of parame- 
ters than in the standard model with three right-handed 
neutrinos. In fact, various recent works in the minimal 
SO(IO) models with a single 10 and a single 126 Higgs 
multiplets0, |M Q> IE 1^ have substantiated this point of 
view. Unlike the generic seesaw models based on the ex- 
tension of the standard model as well as right-handed 
neutrino extended SU(5) model, the above SO(IO) mod- 
els have fewer parameters and lead to predictions for neu- 
trino mixings and phases in their minimal version with- 
out any extra symmetry assumptions. Their predictions 
for solar and atmospheric mixing angles are in agreement 
with present observations Q and that for U^s is not 
far below the present upper limits, making the model 
testable in planned experiments Measurement of Ues 
will therefore provide a crucial test of these models. 

It has further been shown that the simplest way to 
accommodate CKM CP violation in these models is to 
include an additional Higgs field belonging to the 120 
dimensional representation"^. The model with 120 still 



remains predictive in the neutrino sector despite a small 
increase in the number of parameters00 and also leads 
to a solution to the SUSY CP problem. 

Since the quarks and leptons are unified in a GUT, 
it predicts that proton is no longer stable and therefore 
proton decay becomes another test of any GUT such as 
SO(IO). The present experimental bounds on the pro- 
ton lifetime are in fact known to severely constrain some 
GUT models m llH and one must therefore make sure 
that the above class of SO(IO) models are in agreement 
with present experimental results. 

In most generic SUSY GUTs, the dimension five oper- 
ators induced by colored Higgsino,12] provide the domi- 
nant contribution to the proton decay amplitude rather 
than the dimension six operators induced by gauge 
bosons. Since the dimension five operators arise from di- 
agrams involving Yukawa couplings, predictions for pro- 
ton lifetime get related to fermion mass textures. For 
example, in simple GUT theories such as minimal SU(5), 
where by minimal we mean that only the most general 
renormalizable terms are included in the superpotential, 
the up- and down-type quark mass matrices are propor- 
tional to the Yukawa matrices, Yu and Yd respectively. 
In this case the proton decay rate, is directly pro- 
portional to \Yu,ijYd,ki\'^ implying that Tp cannot be arbi- 
trarily adjusted. That is why the minimal SU(5) theory 
is now ruled out by proton decay results [lH. 

One might argue that the minimal SU(5) model is any- 
way not realistic since it predicts wrong relations between 
fermion masses e.g. vris = and — m^. This 
fermion mass problem is however easily remedied in the 
class of minimal SO(IO) modelsQ with a single 10 and 
126 field mentioned above. In view of the fact that it is 
also quite predictive in the neutrino sector 0, 0, 0, 0], it 
is tempting to consider this as the minimal GUT compa- 
rable to what SU(5) GUT was in the 1980's. It is then 
important to look for its predictions for proton decay. 
To be sure, this model like minimal SU(5) (and the stan- 
dard model) does not explain the origin of the Yukawa 
couplings; however if this model is confirmed by experi- 
ments, understanding the Yukawa sector will be the next 
order of business. 
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Unlike the SU(5) model, the A_B = 1 interactions in 
generic SO(IO) modelsHi 

involve more GUT scale sym- 
metry breaking parameters than just the color triplet 
Higgsino mass; therefore the situation for proton decay 
here is less restrictive. In particular, in the minimal 
SO(IO) models of the type discussed in Ref.0], there are 
four free parametersll^ whereas there are about 16 de- 
cay modes which have lower bounds on their partial life 
times. It was shown through a numerical analysis ^3.01 
(without including RRRR operators) that there exists 
a very small region in these parameter spaces for LLLL 
operators, where all the present experimental constraints 
are satisfied for lower tan/3. 

In this paper, we study the proton decay constraints 
on SO(IO) models with 10 + T26 + 120 Higgs fields. 
We show that present proton decay constraints combined 
with fermion masses and mixings imply a specific rela- 
tion among SO(IO) breaking vacuum expectation values 
(VEVs) and a very specific form for the Yukawa tex- 
tures. Roughly, they imply that the proton decay op- 
erators are proportional to the product of two up-type 
quark Yukawa couplings i.e. Yu^ijY^^ki instead of YdYu as 
in SU(5). As a consequence, proton decay is not only 
suppressed for the left-handed quark contributions but 
also for right-handed ones and furthermore, the suppres- 
sion works also for large tan /3 values as well which makes 
this model easily distinguishable from other simple GUT 
models. In addition, it leads to definite predictions for 
the neutrino mixing parameters Ue3 and the Dirac phase 
Smns- We also show the essential role played by the 120 
in the suppression of proton decay when all operators are 
taken into account. These are the new results of this pa- 
per, which have important implications for the viability 
of the minimal SO(IO) model for neutrino masses with 
CKM CP violation. 



FERMION MASS AND PROTON DECAY IN 
MINIMAL SO(IO) 

We start by introducing the Yukawa interactions and 
the contents of Higgs fields in the SO(IO) model. The 
Yukawa superpotential involves the couplings of 16- 
dimensional matter spinors ipi {i denotes a generation 
index) with 10 (H), 126 (A), and 120 (D) dimensional 
Higgs fields: 
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Wy = -h,ji)^ijjjH + -f.jtJjiTpjA + -h[jip,^jD. (1) 

The Yukawa couplings, h and / are symmetric matrices, 
whereas h' is an anti-symmetric matrix due to SO(IO) 
symmetry. One 126 Higgs multiplet A is also introduced 
as a vector-like pair of A whose VEV reduces the rank 
of SO(IO) group. This helps to keep supersymmetry un- 
broken down to the weak scale. In order to break SO(IO) 
symmetry down to the standard model, we employ one 



210 Higgs field (4>)'T6^ which also contains a pair of 
Higgs doublets ^d). Altogether, we have six pairs 
of Higgs doublets: ipd = Z?^, D^Ad, A^, $d), ip^ = 
(i/^",!)^,!?^, A„, A„,$,J, where superscripts 1, 2 of 
stand for SU(4) singlet and adjoint pieces under the 
G422 = SU(4)xSU(2)xSU(2) decomposition. The mass 
term of the Higgs doublets is given as (iPd)n (Mr))nh(v>ii.)h, 
and the expression of the matrix Mo is given in Ref. [ij] . 
The mass matrix of the Higgs doublets is diagonalized 
by unitary matrices U and V: UMdV^ = M'^''^. We 
assume that (M]^"*^)ii = /i, where ^ is a Higgsino mass 
in the MSSM and the mass scale is much smaller than 
the GUT scale. Since we concentrate on the structure of 
Yukawa couplings, we do not need to specify the dynam- 
ical reason of the mass hierarchy in this letter, but we 
just require a fine-tuning such as in the case of minimal 
SU(5) model[i^. The MSSM Higgs doublets are given 
as linear combinations: Hd = U^Jipd)a, Hu = V{^{(pu)a- 
We use "Y diagonal basis" (or SU(5) basis) to de- 
scribe the standard model decomposition of the SO(IO) 
representation^^! ^M- The expression of the Yukawa 
interaction under the 6422 decomposition can be de- 
rived from Ref.0]. The decomposed Yukawa interac- 
tions which give fermion masses are written as 

W^"""^- ^ /ii?]"(g# + ^e") + hHl°{qu'' + iv"") (2) 



+ h'D^diqd" + £e') + h'Dl{qu'' + iv'') 
1 



V3 

where q, u'^, d*^, £, e'^, z^*^ are the quark and lepton fields 
for the standard model, which are all unified into one 
spinor representation of SO (10). We obtain the Yukawa 
coupling matrices for fermions as 



Yu = h + r2f + r3h', 

Yd = n{h + f + h'), 

Ye = riCh-3f + Ceh'), 

Yv = h-3r2f + c„h\ 



(3) 
(4) 
(5) 
(6) 



where the subscripts m, d, e, v denotes for up-type quark, 
down-type quark, charged-lepton, and Dirac neutrino 
Yukawa couplings, respectively and 

h = Vnh, ri = Un/Vn, = nVi^/Ui^, (7) 
iVi2 - Vi3lVi)l{Ui2 + C/i3/\/3), (8) 
/ = 1/14/(73 ri)/, h' = {U12 + Ui3/Vs)/nh', (9) 

Ce = (f/12 - VSUl3)/{Ui2 + C/13/V3), (10) 
iVi2 + \/3Fi3)/(f/l2 + Uls/VS). (11) 

The Majorana mass matrices for both left- and right- 
handed neutrinos are proportional to the coupling /. In 
this letter, we will be using type II seesaw p^. 
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Next we consider dimension five operators induced by 
Higgs triplets. The dimension five operators {LLLL and 
RRRR operators), 

-w,^ \ct'qmmtj + e^«d,^ (12) 

are obtained by integrating out the triplet Higgs 
fields, tpj^ = (7%, D^^, D^, Aj^, AjT, A^, %) and tpr = 

{Ht, Dt, D'rp, At, At, Aj^^^t)- The quantum num- 
bers under SU(3)xSU(2)xU(1)f of the field ipT is 
(3,1,-1/3). In the expression of (pT, the fields with 
"' are decuplet, and the others are sextet under SU(4) 
decomposition. The RRRR operator, Cr, is also gen- 
erated by other colored triplet, ip^ = {D^, A^) and 
ipc = iDc,Ac), where the quantum number of (pc is 
(3,1,-4/3). The mass term of the Higgs triplets are 
given as i(pjT)a{MT)ab{'PT)b + {'Pc)a{Mc)ab{'Pc)b- The 
mass matrices, Mt and Mc, are 7x7 and 2x2 matri- 
ces respectively, and their explicit forms are given in the 
literature^3- The Yukawa couplings which cause proton 
decay are written as 

H^^"P- = hH^ (qi + u^d") + HHt {^qq + e%^) 

+ /AjT {q£ - u'^d'') + /At {^qq - eS") + ^2/5^ e'^w' 

+ V2h'Djr u^d" + V2h'DL qi (13) 

~V2h'DT e^u" + V2h'DT e^u" 
+ 2/Ac d'^e'' + 2h'D-^ u^u" + 2h'Dc d^e". 

The dimension five operators are written by the Yukawa 
couplings h, f and h' as follow: 

C^f"' = chijhki + xifijfki + X2hijfki + xsfijhki 
+Xih[jhki + x<sh[jfu, (14) 
C]^'"' = chijhki + yifijfki + y2hi-jfki + yzfijhki 
+yih[jhki + ysKjfki + yahijh'ki + Vifi^Ki + VsKjKi 
+y9hafjk + yioKih'jk- (15) 

The common coefficient c is given as c = (Afj7^)ii. We 
obtain the other coefficients as the following: 

(xi,a:2,a:3,a;4,a:5) 

= (M5-/,Af5-/,Mf/,V2A/f3\V2A/5-3i), (16) 

(2/1,2/2, 2/3, y4, 2/5, ye, y7,2/8) 
= (A^5-/-v^Afe-/, -AfgV + V2M6i', "^^f/, (17) 

y2A//iV, -V2Af5V + 2Af6V, V2(-Af2V + ^^31'), 

where subscript T of the matrix Mt is omitted. It is 
important to note that j/3 — —x^, namely their signatures 
are opposite between LLLL and RRRR operator. This 
is derived from the fact that Ht and At have opposite 



D-parity. The coefficients j/g, j/io are generated by ipc, 
and 2/9 y 4(M p^)2i, 2/10 = 4(Af^i)ii. 
The Eas. (|14ll7|l are explicit forms of the dimension five 
operators, but we present here more convenient forms to 
describe proton decay. The Higgs triplet mass matrix 
Mt is diagonalized by two unitary matrices, X and Y , 
as XMtV^ = diag(Afi, A.f2, • • • , Mj). Then we obtain 
the useful formula for dimension five operators 

^T' = y. 7^(^-1'^ + ^«4/ + ^/2X,3/l')^J (18) 

x{Yalh + Ya5f)M, 

Cf = -^{Xalh - Xa4f + V2Xa2h%, (19) 

a 

X {Yalh - {Ya5 - V2Yae).f + V2{Ya3 - Ya2)h')u 

+ (2/9,2/10 terms). 

One can make consistency checks to verify the formula 
by considering specific vacua of the theory. For example, 
in the SU(5) limit, when only one of the colored triplets 
is much lighter than the others, i.e. Mi <C Ma [a ^ 1), 
we can obtain the following relations for the diagonal- 
izing matrices from the explicit form of the Higgs mass 
matrices with the SU(5) vacua in the Ref . [itI Ii8| : 

[/ii = Xn, Vii = Fn, Uu = Xu = 0, (20) 

1^15 : ri5 : >l6 = \/3 : 1 : V2, (21) 

Ui2 : Ui3 ■■ Xi2 ■■ Xi3 - Vi2 : ^13 : Y12 : Y13 (22) 

1 : : V2 : V2, (23) 

As a result, r2 ^ 00 with f 0, 

rs = 0, c, = -1 (24) 

for the Yukawa matrices in Eas. H3l6|l and thus, as ex- 
pected, we get the SU(5) relations, Yu = Yj, Yd = Y^ , 
and the dimension five proton decay operators can be 
written in terms of the Yukawa couplings as 

Cf ^ (r,)., {Y^)ki /A/i . (25) 

One can also obtain the flipped-SU(5) liniit^2l'| in a 
similar manner and find the relations, Yu = Yj- , Yd = 
yj, and 

Cf' ~ iYu),.iYd)ki/Mi, Cf ~ (K)y(i"e)/fc/Afi.(26) 

One can also consider other typical vacua such as the 
one corresponding to C?422- On this vacuum, wc have 
Viz — Vi3 — Ui4 = Vi5 = and all four fermion Yukawa 
couplings are unified as Yu^dje,^ = h + h' , which is not 
viable for phenomenology. In any case, since the doublet 
and triplet Higgs fields are not unified in one multiplet on 
this vacuum, proton decay operator cannot be written in 
terms of fermion Yukawa couplings contrary to the SU(5) 
and flipped-SU(5) limit. The operators C];^^ in the G422 
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limit are symmetric under the interchange of indices ij 
and kl in addition to the symmetry under the exchange 
of individual indices i and j. One can write down the 
dimension five operators as 

CL'^{h + zf){h + zf), Cn^{h-zf){h-zf), (27) 

with z — Xii/Xii. It is important to stress that there 
are relatively opposite signatures for the / contributions 
between Cl and Cr. This is due to D-parity, which we 
have already mentioned. 

SUPPRESSION OF PROTON DECAY 

Let us now investigate the conditions required to sup- 
press the proton decay rate in the minimal SO(IO) model 
described above. We first note that the four-fermion pro- 
ton decay operators are produced by gaugino and Hig- 
gsino dressing of the dimension five operators. The four- 
fermion operators of LLLL type get dominant contribu- 
tion from wino exchange and therefore retain the same 
flavor structure as that of the original dimension five su- 
persymmetric operator (we have assumed universality of 
scalar masses for suppressing FCNC processes). How- 
ever, as far as the RRRR operators are concerned, they 
receive contributions only from Higgsino exchange which 
can involve the top quark and r lepton Yukawa couplings. 
For instance in the case of SU(5) model, the contribution 
of the RRRR operator to the decay mode p Kvr gets 
so enhanced because of this that it exceeds the current 
experimental bound for any tan (3 as long as stop mass is 
less than around 1 TeV[ll|. 

Since the original Yukawa couplings h, /, and h' are 
functions of fermion Yukawa couplings 1^, Yd and Yg via 
Eqs.lPElj their textures are roughly determined from the 
experimental inputs of quark and lepton masses and mix- 
ings. For instance, to fit the strange quark and muon 
masses, bottom-tau unification, relations among CKM 
mixings and also large mixings of neutrinos, the cou- 
pling / is almost determined to have the form / ^ 
/A2 A A\ 

All ms/mii, where A ~ 0.2. A naive impli- 
V A 1 1/ 

cation of this is that since up-type quark masses are 
more hierarchical than down-type ones (i.e. niu/mt ^ 
md/rub, rric/rnt ^ ms/m},), the expression for the up- 
type Yukawa matrix, Y^, in Eq.Q requires the following 
two typical choices: (a) there is cancellation among /i, /, 
and h' , or (b) h itself has a hierarchical form similary to 
the up-type quark masses. The first choice (a) is the case 
where [1,2] block of hij is not far smaller than /y, but 
^■2,3 are chosen to be certain values to make rrin, nic are 
hierarchically small. The second choice (b) corresponds 
to the case where r2,3 ^ 0. The second choice appears to 
be required to suppress the proton decay. Let us discuss 
the reason. 



In order to suppress the decay rate, we need small cou- 
plings for first and second generations in the expressions 
in Eas. H14ll5(l . Clearly this would also require a cancel- 
lation among h, f and h' if we take the first choice (a). 
Since in general the coefficients r2 , rs in up- type Yukawa 
matrix and Xi and yi in proton decay operators are unre- 
lated, one must find a situation where both cancellations 
can be achieved in a satisfactory manner so as to be con- 
sistent with all data. However, if we take the choice (a), 
the cancellation cannot happen naturally due to the fol- 
lowing reasons. 

First, let us discuss the 126 Higgs contribution. Since 
there is an opposite signature for one of the coefficients in 
the LLLL and RRRR operators, ys — —X3, the cancella- 
tion required to obtain small Yukawa coupling for Yu by 
tuning r2f cannot simultaneously suppress both LLLL 
and RRRR operators. Thus, in general, it is hard to 
suppress proton decay rate by tuning 126 colored Higgs 
mixing X14. Next let us see the contribution from 120 
Higgs field. The coefficient Cf^''^ is symmetric in the in- 
dices kl due to SU(3)xSU(2) contraction. Therefore h' 
contribution is absent for the qkqi part in LLLL opera- 
tor, whereas it is of course present in the fermion masses. 
Thus, if the cancellation in Yu happens by tuning r^h' , 
such cancellation will not help in suppressing the LLLL 
operator. 

The above discussions lead to the fact that the proton 
decay rate cannot be suppressed in natural way if we take 
the choice (a). As we will see numerical studies later, we 
require a fine-tuning to the level of 0.01% to suppress all 
the proton and neutron decay modes when we consider 
general parameter fitting in the choice (a) . The choice (b) 
(r2,3 ~ 0) is necessary to achieve natural suppression of 
proton decay as a result. 

We now show how the proton decay rate in the 
choice (b) is suppressed rather than the minimal SU(5). 
If r2,3 ~ 0, the LLLL dimension five operator can be 
a form as Cf^'^'' cx {Yu + ^h')ij{Yu)ki and in the RRRR 
operator C^^^\ kl part is also related to Yu- This will 
correspond to the case where Xi4,Yi5 ^ 0. We will 
give an example later where this can happen but first 
we note that for the case just noted, the RRRR con- 
tribution to p —i Kvr mode is suppressed compared to 
the usual SU(5) models for the entire range of tan/3 up to 
50. Here, the dominant contribution is proportional to Au 
giving a suppression factor Au/A^ ^ 1/100 for tan/3 ~ 50 
compared to minimal SU(5) case. Similarly, since the kl 
part of Cl are also related to the Yu instead of Y^, the 
LLLL contribution to the p — > KD is also suppressed 
even for tan j3 ^ 50, compared to the SU(5) model (since 
Ac/As ~ 1/5)- However as it turns out, these suppres- 
sions are not enough. We need to specify the Yukawa 
texture for the purpose as we see below. 

Before describing the specific choice of Yukawa tex- 
tures, let us show the numerical values for the pro- 
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TABLE I: Typical proton decay amplitude for different values 
of r-2 and r^. To satisfy current experimental bound naturally, 
\Ahh\ < 10"* is needed. 



r2, r-3 


Ahh 


r2 = r3= 
^■2= 0, 7-3= 1 
r2= 0.5, r3= 


(4.2 + 0.083i) X 10^^ 
(1.4 - 0.17i) X 10"" 
(7.3 + 2.6i) X 10"^ 



ton decay amplitude to confirm that r2,3 ~ is nec- 
essary to satisfy the current experimental bounds nat- 
urally. Rewriting the proton decay amplitude as A — 
o-iPp/ {4^TTMTmsusY)A, we define 

A = cAhh+xiAff+x2Ahf+x3Afh + XiAh'h + - ■ ■ ■ (28) 

The coefficients c and Xi are given in Ea. ll4|l . The 
RRRR contribution can be written in the same way. To 
satisfy the current proton and neutron decay bounds, 
we need \Ap^Ki^\ S 10"*, |A„^^p| < 2 • 10"* and 
\An^Kv\ ^ 5 • 10"* if the colored Higgsino mass is 2 • 10^^ 
GeV, and squark and wino masses are around 1 TeV and 
250 GeV, respectively. In order to satisfy those bounds 
naturally, we need Ahh ^ 5 • 10~*. If Ahh ^ 10"^, we 
need to tune Xi and yi for every decay mode to can- 
cel Ahh, which is unnatural. (Further, assuming c — > 
cannot make successful suppression of the decay ampli- 
tude since in that assumption we also need to suppress 
Aff which can only be suppressed for tan/3 < 3). The 
Ahh depends on the magnitudes of the elements from the 
[1,2] block of h which is determined from the fit to the 
up- type Yukawa as a function of r2 and r^,. In Tabled we 
show the Ahh component in the proton decay amplitude 
for different values of and r2- The Ahh is calculated 
in the basis where Yu — diag(AM, Ac, At). As we can see 
from the table, that Ahh can be very large (for large r2,3) 
and this requires very high level of fine-tuning for all the 
decay modes (including charged leptonic decay) to sup- 
press nucleon decay as we noted. The smallness of r2 
and ra is required to suppress Ahh- However, even in 
the case r2 = — 0, the Ahh amplitude is not enough 
suppressed {Ahh ~ 10"^) and we still need fine-tuning 
among the coefficients to satisfy the data. This is mainly 
because A„_c are not enough small to satisfy Ahh ^ 10"*. 

We now present the specific choice of Yukawa texture 
to suppress the proton decay rate further. The texture 
is given as h ~ diag(0, 0, 0(1)), so that up and charm 
quark masses arise from contributions of / and h' : 

/ ^ ~ A3 \ / A3 A3 \ 

f^\^Q\^X^\,h'^i -A3 AM, 
V A3 A2 A2 / \ -A3 -A2 / 

(29) 

where A ~ 0.2. In this texture, = 0, and r2 is fixed as 
Tivnajm^i ~ Ac {r-i — 0.1) to generate the correct charm 



mass. We have /11.12 ^ (nonzero value of r2 is there- 
fore allowed) to ensure small values of hw^yi, and h'yi 
generates the down-quark mass and Cabibbo angle 
with md/ms — sin^ 9^- The up-quark Yukawa coupling 
is generated as A^, ~ {r2^^Y , and we also have a rela- 
tion mdmsmh ~ clmf.m^mr, where ~ 1 in the pre- 
ferred vacuum for the 120 Higgs coupling, Ea. (|24|l . In 
the basis where is diagonal, Ahh in this texture is not 
completely zero but can become much smaller than 10"*. 

We comment that we also need to examine the contri- 
bution of the other components e.g. Af f^hfjh.h' f.h'h.--- 
after we suppress the Ahh- However, their contribu- 
tions are associated with the colored Higgs mixing angles, 
which are inputs and can be suppressed by our choice of 
the vacuum expectation values and the Higgs couplings. 
According to our numerical studies, some of the mix- 
ing angles must be about a few percent in the case of 
tan/3 ~ 50 to suppress the decay. However, the mix- 
ing angles can become larger as tan/3 becomes smaller. 
Details will be given in a future paper [2^. 

Without the particular choice of texture, as mentioned 
above, the proton decay can not be suppressed naturally 
in these models unless the tan /3 is very small. 

The proton life time for p — > KD for this choice of tex- 
ture can be larger than the current experimental bound, 
Tp > 2 • 1033 years for any tan/3 (using the lightest col- 
ored Higgsino mass to be 2 • 10^^ GeV and squark mass 
scale around 1 TeV). All other nucleon decay modes are 
suppressed as well. In our calculation, we use long- and 
short-distance renormalization factor, Al = 1.43 and 
typically ^5 = 1.8 similar to Ref.ji^. 

Given the above texture for Yukawa couplings, we find 
that Ue3 is restricted to a range. In figure ^ we plot 
Ue3 as a function of r2 with /ii,i2 — * (as required in 
the most preferred texture) and the value at \r2 \ ~ 0.1 is 
most important. We find that Ues is between 0.07 to 0.09. 
For this fit, sin^ 26*23 is maximal and tan^ 6'soiar ~ 0.4. 
The prediction for sin|(5A/Ars| lies between 0.3 to 0.7 for 
\r2 \ — 0.1 as shown in figure El The Yukawa matrices are 
assumed to be hermitian in order to keep the model free 
from SUSY CP problem. 

The presence of h' is a necessity to suppress proton 
decay (suppress Ahh) and fit the fermion masses. This 
h' also helps to explain CKM CP violation^. 

Finally, we show how the above proton decay suppres- 
sion arises by an adjustment among different VEVs (or 
symmetry breaking parameters). In the above, we need 
to have ra ~ in the Eq.®. Since this is satisfied in 
Ea. (|24|l for the SU(5) condition, it may be a hint that 
we stay close to the SU(5) symmetric vacuum. Secondly, 
we need r2 — while suppressing 126 colored Higgs con- 
tribution to proton decay, namely only Uu is enhanced 
while other mixings for the sub-multiplets in 126 are 
small, 

Uu » Vi5, Xi4, Yi5- (30) 
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where 77 is a coupling of $AA term. We also have similar 
terms for V15, Y15 and Yie- AH these terms have different 
denominators. All the Higgs mixing angles tend to zero 
in the limit A2 0. However, suppose that •\/65'i + 
\/2$2 + 8<i>3 ~ is satisfied, only Uu can be of finite 
value, all other mixing angles are zero and Ea. l|3U|) is 
satisfied. On such a vacuum, the proper strange quark 
and the muon masses can be realized by an enhancement 
of [/i4, and Yukawa coupling of up- type quark, Yu, is 
almost proportional to h. This is just an example, and in 
our detailed quantitative work, we keep all other terms in 
the Higgs potential and we satisfy Eas. H24l30|) to suppress 
the proton decay rate. 



''2 



FIG. 1: \Ue3\ is plotted as a function of r2- 
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In conclusion, we have analyzed the fermion masses 
and dimension five AB = 1 operators in the minimal 
SO(IO) model with 10, T26 and 120 Higgs fields cou- 
pling to matter. We show that by a choice of suitable 
textures, one can not only get correct fermion masses 
and mixings but also suppress the contributions to pro- 
ton decay from both the LLLL and RRRR operators 
for the entire range allowed tan/3 parameter of MSSM. 
This choice of textures requires a suitable SO(IO) break- 
ing vacuum condition which is close to SU(5) invariant 
vacua. In the most favorable region of parameter space 
Ue3 is predicted to be 0.07 to 0.09 and sinlJ^/ArsI to be 
0.3 to 0.7, which can be used to test the model. 



FIG. 2: sin \5k 



is plotted as a function of r2- 
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In order to explain how these conditions may be satisfied 
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the following relation: ^/Q<^i = \/2$2 = ^3 (We have 
used the same normalization as in the Ref.[l7|). They 
lead to the SU(5) relations H20I23|I . Now we perturb the 
Higgs potential with a small couphng, A2-ff A$. We ob- 
tain r3 oc Ai(-\/6<I'i — *I'3) (where Ai is associated with 
XiHD^ term). If 76$! = $3, we have rs ~ 0. We 
obtain the Higgs mixings 



C/i4Ci -6%/5A2/ry 



^/2^2 - *3 



Xi4~ -2yT5A2/?7 



Ve^i + \/2$2 + 8$3 
%/6$i - %/2$2 



V6$i + 3\/2$2 + 6$3 
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